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Abstract 

We consider a class of nonlinear differential equations that arises in the 
study of chemical reaction systems that are known to be locally asymp- 
totically stable and prove that they are in fact globally asymptotically sta- 
ble. More specifically, we will consider chemical reaction systems that are 
weakly reversible, have a deficiency of zero, and are equipped with mass 
action kinetics. We show that if for each c G W^q the intersection of the sto- 
ichiometric compatibility class c + S with the subsets on the boundary that 
could potentially contain equilibria, L\y, are at most discrete, then global 
asymptotic stability follows. Previous global stability results for the sys- 
tems considered in this paper required {c + S) H Lyy = for each c G IR>0' 
and so this paper can be viewed as an extension of those works. 

1 Introduction 

This paper is motivated by the consideration of a class of nonlinear systems that 
arises in the study of chemistry and biochemistry. Suppose there are m chemical 
species, {Xi, . . . , Xm}, undergoing a series of chemical reactions. For a given re- 
action, denote hy y,y' E Z>q the vectors representing the number of molecules of 
each species consumed and created in one instance of that reaction, respectively. 
Using a slight abuse of notation, we associate each such y (and y') with a linear 
combination of the species in which the coefficient of Xi is yi. For example, if 
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y = [1, 2, 3]^ for a system consisting of three species, we associate with y the 
linear combination Xi + 2X2 + 8X3. Under this association, each y (and y') is 
termed a complex of the system. We may now denote any reaction by the notation 
y y', where y is the source, or reactant, complex and y' is the product complex. 
We note that each complex will typically appear as both a source complex and a 
product complex in the system. Let S = {Xi},C = {y}, and IZ = {y ^ y'} 
denote the sets of species, complexes, and reactions, respectively. Denote the 
concentration vector of the species as x G M™. In order to know how the state of 
the system is changing, we need to know the rate at which each reaction is tak- 
ing place. Therefore, for each reaction, y y', there is a function Ry^yi{-) 
satisfying the following: 

1. Ry-,yi{-) is a function of the concentrations of those species contained in 
the source complex, y. 

2. Ry^yi{-) is monotone increasing in each of its inputs and Ry^yi{x) = if 
any of its inputs are zero. 

The dynamics of the system are then given by 



where the last equality is a definition. The functions Ry^y' are typically referred to 
as the kinetics of the system. This notation closely matches that of Feinberg, Horn, 
and Jackson, and it is their results that the main theorem in this paper extends 



Therefore, x(t) — x{0) remains in the linear space S = Span{y' — y}y^y'en for 
time. We shall refer to the space S as the stoichiometric subspace of the system 
and refer to the sets c+ S, for c G W^, as stoichiometric compatibility classes, or 
just compatibility classes. Later we will demonstrate that trajectories with positive 
initial conditions remain in W^q for all time. The sets {c + 3)0 W^q will therefore 
be referred to as the positive stoichiometric compatibility classes. Given that tra- 
jectories remain in their positive stoichiometric compatibility classes for all time. 




(1) 



(ffmiTiiiiiioiii). 

Integrating equation ([T]) gives 



x{t)=x{0)+ J2 ( f Ry^y'i<s))ds)iy'-y) 



2 



we see that the types of questions that one should ask about these systems dif- 
fer from the questions one normally asks about nonlinear systems. For example, 
instead of asking whether there is a unique equilibrium value to the system Qi, 
and then asking about its stability properties, it is clearly more appropriate to ask 
whether there is a unique equilibrium within each positive stoichiometric compat- 
ibility class and, if so, what are its stability properties relative to its compatibility 
class. 

The most common kinetics chosen is that of mass action kinetics. A chemical 
reaction system is said to have mass action kinetics if 

for some constant ky^y/. It has been shown that for many systems of the form ([T]) 
with mass action kinetics, there is within each positive stoichiometric compatibil- 
ity class precisely one equilibrium and that equilibrium is locally asymptotically 
stable relative to its class ( [fTTl lQllSlI'). In order to show that the equilibrium values 
are locally stable, the following Lyapunov function is used (one for each compat- 
ibility class) 

m 

V{x, x) = V{x) = ^ [xi(\n{xi) - \n{xi) - 1) + Xi] , (3) 

i=l 

where x is the unique equilibrium of a given positive stoichiometric compatibility 
class. It turns out that the function V "almost" acts as a global Lyapunov func- 
tion. That is, V is non-negative for x E {x + S)ri M™q, zero only at x, and strictly 
decreases along trajectories. However, V does not tend to infinity as trajectories 
near the boundary of (x + S) PiR^q, and without such unboundedness one can not, 
in general, conclude global stability. It has been shown in numerous papers, how- 
ever, that global stability of x does hold if there are no equilibria on the boundary 
of (x + S) nMII'o (IH mini ID- Therefore, work has been done giving sufficient 
conditions for the non-existence of boundary equilibria in order to conclude that 
the equilibrium value within each compatibility class is globally stable relative to 
its class (E mill). 

To each subset W of the set of species, the set of boundary points is 
defined to be 

Lw = {x eW" : Xi = ^ Xi e W}. (4) 
We will show that there are no boundary equilibria if and only if 

[(c + ^) n M^o] n = 0, (5) 



3 



for all c G M>Q and for certain subsets of the species, W. We will then prove that 
global stability holds if the intersection given in equation © is either empty or 
discrete for each c G ]R™o those same subsets, W. This will imply that global 
stability holds even if there are boundary equilibria, so long as the boundary equi- 
libria are extreme points of the positive stoichiometric compatibility classes. To 
the best of our knowledge there is only one other result concerning the global sta- 
bility of mass action systems with boundary equilibria, and is contained within the 
Ph.D. thesis of Madalena Chavez ([6J). In order to guarantee global stability even 
if there exist boundary equilibria, she requires that each boundary equilibrium be 
hyperbolic with respect to its stoichiometric compatibility class and she requires 
another (more technical) condition on the stable subspaces of each boundary equi- 
librium (see 1161, pg. 106 for details). As our result are applicable to systems with 
boundary equilibria that are non-hyperbolic, our results can be viewed as an ex- 
tension of those in [jH. 

The layout of the paper is as follows. In Section [21 we will introduce the 
systems we consider in this paper: weakly reversible, deficiency zero systems 
with mass action kinetics. We will then present some preliminary results and will 
conclude with a proof that global stability follows if there are no equilibria on 
the boundary of the positive stoichiometric compatibility classes. No originality 
is claimed for this result as it is known. Also in Section l2l we demonstrate how 
the "no boundary equilibria" assumption is equivalent to equation © holding for 
all c G M™Q and certain subsets of the species, W. In Section l3l we extend the 
previous theorems to prove that global stability still holds if the intersection given 
in equation Q is always either empty or discrete for those same subsets, W. We 
also show in Section l3] how the hypothesis that the intersection in equation Q is 
always empty or discrete is equivalent to a condition on the extreme points of the 
non-negative stoichiometric compatibility classes. In Section IH we demonstrate 
our results on a number of examples. Finally, in Sectionl^we sketch how to extend 
our results to systems with non-mass action kinetics. 

2 Preliminary results 

We start with definitions taken from lUD, lH, and lH. 

Definition 2.1. A chemical reaction network, {S, C, TZ}, is called weakly re- 
versible if for any reaction y y', there is a sequence of directed reactions be- 
ginning with y' and ending with y. That is, there exists complexes yi, . . . ,yk such 
that the following reactions are inTZ: y' ^ yi, yi y2, yk ^ V- 
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To each reaction system, {S,C, TZ}, there is a unique, directed graph con- 
structed in the following manner. The nodes of the graph are the complexes, 
C. A directed edge is then placed from complex y to complex y' if and only if 
y ^ y' E TZ. Each connected component of the resulting graph is termed a 
linkage class of the graph. We denote the number of linkage classes by /. 

As shown in the introduction, each trajectory remains in its stoichiometric 
compatibility class for all time. There is another restriction on the trajectories of 
solutions to equation ([T]) that is given in the following lemma. The proof can be 
found in both [14.1 and UJ • 

Lemma 2.2. Let x{t) be a solution to ^ with initial condition a;(0) G M^q. Then, 
x{t) G W^^Jorallt > 0. 

2.1 Persistence and cj-limit points 

By Lemma [Z2l each trajectory must remain within W^q if its initial condition is 
in ffi>0' therefore the linear subsets of interest are the intersections of the stoi- 
chiometric compatibility classes and M^q. Recall that in the Introduction these 
sets were termed the positive stoichiometric compatibility classes. This paper will 
mainly be concerned with showing that trajectories to systems given by equation 
([B remain away from the boundary of the positive stoichiometric compatibility 
classes. That is, we will show that the systems are persistent. To be precise, let 
(j){t,^) be the solution to equation ([T]) with initial condition ^ G ]R™q. The set of 
a;-limit points of the trajectory is 

"^(0 = E M>o : (f){tn,0 — > X, for some tn — >■ oo}. (6) 

Definition 2.3. A system is persistent if uj{^) fl dW^Q = 0, for each ^ G M>o- 

We refer the reader to f'?! 121 [H [151 for some of the history and usage of the 
notion of persistence in the study of dynamical systems. In order to show that 
a chemical system is persistent, it is critical to understand which points on the 
boundary are capable of being tu-limit points. With that in mind, we introduce the 
following definition. 

Definition 2.4. A nonempty subset W of the set of species is called a semi-locking 
set if for each reaction in which there is an element of W in the product complex, 
there is an element of W in the reactant complex. W is called a locking set if 
every reactant complex contains an element of W. 
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Locking and semi-locking sets are easily understood. Suppose that W C 
{Xi, . . . , Xm} is a locking set. Then, because every reactant complex contains 
an element of W, if the concentration of each element of is zero, each ki- 
netic function, Ry~,y', must equal zero. Therefore, all of the fluxes are zero, and 
x{t) = 0. We therefore see that the system is "locked" in place. Now suppose 

is a semi-locking set. If the concentration of each element of W is zero, then 
any flux which affects the species of W is turned off and the elements of W 
are "locked" at zero. Semi-locking sets have another, important, interpretation in 
terms of the linkage classes and weak reversibility. If the concentrations of the el- 
ements of a semi-locking set are equal to zero and the system is weakly reversible, 
then all of the fluxes of any linkage class with a complex containing an element of 
W are equal to zero (and so these linkage classes are "locked"), while the fluxes 
of the other linkage classes are all not equal to zero. Therefore, the concept of 
a semi-locking set and a locking set are equivalent for systems that are weakly 
reversible and have only one linkage class. We note that our notion of a semi- 
locking set is analogous to the concept of a siphon in the theory of Petri nets. See 
O for a full discussion, including historical references, of the role of Petri nets in 
the study of chemical reaction networks. 

The following Theorem now characterizes the boundary points that have the 
capability of being cu-limit points of the system. This result was first proved in 
[|2l, however the proof given here is completely different and straightforward. 

Theorem 2.5. Let W be a non-empty subset of the species. If there exists a C, E 
M™Q such that uj{^) n Lw 7^ 0, then W is a semi-locking set. 

Proof. Suppose, in order to find a contradiction, that there is a ,^ G M^g ^^'^ ^ 
subset of the species, W, such that u;(^) fl Lw 7^ and W is not a semi-locking 
set. Lety G uj{^)r\Liv. We note that there exists a species , withX^ G W, such 
that at least one input to Xj (term in fj of equation ([T]) with a positive coefficient) 
is non-zero if the concentrations are given by y, for otherwise W would be a semi- 
locking set. Therefore, because all outputs from species Xj (terms in fj with a 
negative coefficient) are zero at y, there exists e > and k > such that if 
x{t) G M™o n B,{y), then 

f,ix{t)) = x'^it) > k, (7) 

where B^(y) = {x : \x — y\ < e}. 

Because /(■) is C\ we have ||/||oo,«oc < M for some M > 0, and this bound 
is valid in M^^ n B,{y). Therefore, for any < a < 6, if x(t) G M™o ^ B,{y) for 
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t G (a, b), we have that 



\x{b) — x{a)\ 



b 



f{x{s))ds 



<{h- a)M. (8) 



Now consider a partial trajectory starting on the boundary of R™q fl B^{y) at time 
te, ending on the boundary of M^q fl B^/2{y) at time and remaining within 
that annulus for all time in (te,te/2)- Note that one such partial trajectory must 
exist every time we enter M™q fl 5^/2 (y), and this happens at least once by our 
assumption that y e a;(^) fl Lw- By dS]), t,/2 -te> e/(2M). On the other hand, 
by ©, x'j{t) > kfoxte (te, te/2)- Therefore, 

Xj{te/2) = XjiQ + x'As)ds 



> Xj{Q + ek/{2M) 

> ek/{2M). 

Combining the above with the fact that we still have x'j{t) > k on ]R™q fl B^/2{y), 
and we see that there can not exist times t„ such that ^ y, as n ^ oo. This 
is a contradiction and completes the proof. □ 

Remark. Theorem [23] is a powerful tool for understanding the dynamics of chem- 
ical reaction systems. We see that in order to prove that a chemical system is 
persistent, it is sufficient to show that [(c + 5) n M>o] n Lw = ^ for all c G M>o 
and all semi-locking sets W. We will show in Lemma [2^ that for many reaction 
systems such a condition is equivalent to having no equilibria on the boundaries 
of the positive stoichiometric compatibility classes. 



2.2 Deficiency and the Deficiency Zero Theorem 

We require one more definition before we can state precisely the types of systems 
we consider in this paper. 

Definition 2.6. The deficiency, 5, of a reaction network {S, C, TZ} is given by 
5 = n — I — s, where n is the number of complexes of the system, / is the number 
of linkage classes, and s =dim S, the dimension of the stoichiometric subspace. 

Remark. It has been shown that the deficiency of a reaction network is a non- 
negative number. In fact, the deficiency is the dimension of a certain subspace 
associated with the system. See [HI, [[TOl , and [[U for details. 
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The main types of systems considered in this paper are those with mass action 
kinetics, are weakly reversible, and have a deficiency of zero. The following 
theorem by Feinberg (|[8l[T0l) is the catalyst for studying such systems. The proof 
can be found in [3 or [fTOl . 

Theorem 2.7 (The Deficiency Zero Theorem). Consider a system of the form ([T]) 
with mass action kinetics that is weakly reversible and has a deficiency of zero. 
Then, within each positive stoichiometric compatibility class there is precisely one 
equilibrium value and it is locally asymptotically stable relative to its compatibil- 
ity class. 

In order to prove that the systems considered in the Deficiency Zero Theorem 
are locally asymptotically stable, the Lyapunov function ([3]) is used. It is shown 
that for X G (x + 5) n M!!?o (where x is the equilibrium guaranteed to exist by The- 
orem l2.7l) . V{x) > with equality if and only if x = x, and that dV{x{t))/ dt < 
for all trajectories with initial condition in {x + S) H M!!?q. We will make use of 
these facts throughout the paper without reference, however we point the inter- 
ested reader to the original works for details. 

2.3 Boundary equilibria 

The following Lemma shows that having no equilibria on the boundaries of the 
positive stoichiometric compatibility classes is equivalent to [(c + S") fl M>q] fl 
L\y = for all c G M™q and semi-locking sets W. Following Lemma [Z8] we 
present a theorem pertaining to any system equipped with a globally defined Lya- 
punov function that does not necessarily go to infinity as x goes to the boundary 
of the domain. We then use these results in combination with Theorem 12.51 to 
conclude that for weakly reversible, deficiency zero systems with mass action ki- 
netics, having no equilibria on the boundary of the positive stoichiometric compat- 
ibility classes implies global stability of the equilibrium values given by Theorem 
12.71 We again note that it is already known that global asymptotic stability fol- 
lows from a lack of boundary equilibria. For example, in [14] Sontag showed that 
all trajectories must converge to the set of equilibria, and so a lack of boundary 
equilibria implies convergence to the unique equilibrium in the interior of the pos- 
itive stoichiometric compatibility class. We rederive this result here because our 
methods put it in a larger context in which global stability is understood through 
the intersections given in equation ^ and because it makes clear how our results 
in Section |3] are truly a generalization of this fact. 
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Lemma 2.8. For any chemical reaction system, the set of equilibria on the bound- 
aries of the positive stoichiometric compatibility classes is contained in |J^y [(c+ 
S) n W^q] n Lw, where the first union is over c € M™g and the second union is 
over the semi-locking sets. Further, if there are no equilibria on the boundaries 
of the positive stoichiometric compatibility classes for a weakly reversible, defi- 
ciency zero system with mass action kinetics, then [(c + S*) fl M>o] ^ = ^for 
all c G M™Q and semi-locking sets W. 

Proof. Let y be an equilibrium on the boundary of a positive stoichiometric com- 
patibility class. Let W be the set of species with a concentration of zero at y. 
Because each complex that contains an element of W is providing zero flux, in 
order for y to be an equilibrium value each reaction in which there is an element 
of W in the product complex must have an element of W in the reactant complex. 
Thus, W is a semi-locking set and y G [(c + S") Pi M^q] ^w-^ for some c G K>0' 
In order to prove the second part of the Lemma, we suppose is a semi- 
locking set for the system and suppose y G [{c-\-S) nR>o] I~'-^vk for some c G M™o- 
We will now produce an equilibrium value on the boundary of {c + S) Ci M.^q. If 
W = S, y = 0, and, because is a semi-locking set, y is an equilibrium. 
Otherwise, consider the system consisting only of those species not in the semi- 
locking set W. By the arguments in ifTOl . the linkage classes not "locked" by W 
form their own weakly reversible, deficiency zero system. Therefore, there is an 
equilibrium for that reduced system, z. Let y = (z, 0) (where we have potentially 
rearranged the ordering of the species so that those not in the semi-locking set 
came first), y is a boundary equilibrium value to our original system. Therefore, 
the result is shown. □ 

Theorem 2.9. Let x{t) = x{t,x{0)) be the solution to x = f{x) with initial 
condition x{0), where f is and the domain of definition of the system is the open 
set C C M™. Let x E C be the unique equilibrium value to the system. Finally, 
suppose that there is a globally defined Lyapunov function V that satisfies: 

L V{x) > with equality if and only ifx = x. 

2. dV{x{t))/ dt < with equality if and only ifx{t) = x. 

3. V{x) ^ oo, as \x\ oo. 
Then either x{t) x or x{t) — > dC. 
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Proof. Suppose that x{t) x. Because V{-) decreases along trajectories, the 
value V{x{t)) is bounded above by V{x{0)) for all t > 0. Therefore, because 
V{x) ^ oo as |a;| oo, x{t) remains bounded for all t > 0. Also, the local 
asymptotic stability of x combined with the fact that x{t) ^ x implies there is a 
p > such that \x{t) — x| > p for all t > 0. 

Let e > and for x G C let d{x, dC) represent the distance from x to the 
boundary of C. Let = {s G C | d{x,dC) > e and |a; — x| > p}. Using 
that trajectories remain bounded for all time, we may use the continuity of the 
functions W and / to conclude that there is a positive number t] = 77(e) such that 
W{x) - fix) < —T] for all x E C^- Therefore, the maximum amount of time that 
any trajectory can spend in the set is V{x{Qi))/ri (for, otherwise, x{t) x). 
Because e > was arbitrary we see that x{t) dC. □ 

Corollary 2.10. If there are no equilibria on the boundaries of the positive sto- 
ichiometric compatibility classes for a weakly reversible deficiency zero system 
with mass action kinetics, then the unique positive equilibrium value within each 
positive stoichiometric compatibility class is globally asymptotically stable rela- 
tive to its compatibility class. 

Proof. This is a direct result of Theorems 12.51 and 12.91 and Lemma IZ8l □ 

3 Main results 

By Lemma 12. 8[ we see that the no boundary equilibria assumption for weakly 
reversible deficiency zero systems with mass action kinetics is equivalent to the 
assumption that [(c + 5*) n M>o] n L^/ = 0, for all c G and all semi-locking 
sets W . This then implies global stability by Corollary 12. 10[ We will extend these 
results by proving that global stability holds if [(c + S*) fl ]R>q] fl Lw is empty or 
discrete for each c G M^g Q^ch. semi-locking set W . The following definition 
is necessary. 

Definition 3.1. For a vector x G the support of x, denoted supp(a;), is the 
subset of the species such that Xj G supp(x) if and only if Xi 7^ 0. 

Proposition 3.2. Let {S, C, 71} be a weakly reversible, deficiency zero, mass ac- 
tion chemical reaction system with dynamics given by equation ([T]). Suppose that 
y G uj{x{0)) for some x(0) G M™o- Then there must exist a nonzero zq E S with 
supp{zo) C supp{y). 
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Proof. If y G IR>0' there is nothing to show. Therefore, assume that y is on 
the boundary of the positive stoichiometric compatibility class. By Theorem |2.5[ 
there is a semi-locking set W such that y e [{x{0) + S) n ]R>o] ^ ^w- 
Let V{x) : ]R™q M, be given by equation (|3]) and let 

Vi{xi) = Xi{ln{xi) - In(xi) - 1) + Xi, 

so that V{x) = YllLi Vii^i)- Reordering the species if necessary, we suppose 
W = {Xi, . . . , Xd}. Choose p > so small that for each i < d, Xi < p =^ 
ln{xi) — ln{xi) < 0. Let e > satisfy e < p. Let be a time such that Xi{t^) < e 
for all i < d and \xj{te) — yj\ < e for all j > d + 1. Let = min{t > : 
\xi{t) — yi\ < Xi{te)/2, for alH < m}. We know such and exist because 
y is an u;-limit point of the system. Note that > and that for each i < d, 
Xi{T^) < Xi{te). We consider how V{x{t)) changes from time to time T^. 
Applying the mean value theorem to each Vi(-) term gives 

m 

V{x{T,)) - V{x{Q) = ^ V,{x,{T,)) - V,{x,{Q) (9) 

i=l 
d 

= ^(In(xi) - ln(xi))(xi(T,) - 

(10) 

+ (In(xi) - ln(xi))(xi(T,) - Xi{t,)) 

i=d+l 

for some Xi G [xi{Tt,), Xi{t^)]. Recalling that V decreases along trajectories of 
x{t) by Theorem 12. 7 [ we have V{x{T^)) — V{x(t^)) < 0. Note that because 
for j > (i + 1, we have \xj — yj] < e, there are positive constants Cj such that 
Cj > I In(xj) — Inxj |, and that bound is valid /or any e < p. Let C = J2f=d+i ^j- 
By our choices above, we know that for each i E {1, . . . ,d} the following 
inequalities hold: 

1. In(xi) - In(xi) < 0. 

2. x,(T,) - XiiQ < 0. 

Therefore, each piece of the first sum in equation (flOl ) is strictly positive. Thus, to 
ensure that V is decreasing along this trajectory, the second sum in equation (fTOl) 
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must be negative and, letting Axj = Xi(T^) — Xi(t^) for each i, we have 



In(xi)) Axi < 



i=l 



j=d+i 

m 

j=d+i 



Xj\. 



In(xj)) Ax J 



(11) 



In fact, because each term on the left hand side of equation (fTTI) is positive, a 
similar inequality must hold for each i = 1, . . . , d. That is, fori<d 



(In(xi) - \n{xi))Axi < ^ 

j=d+i 



Xj\. 



For each i < d, Xi E [xi{T^), Xi{t^)] and Xi{T^),Xi{te) < e. Hence, letting 
I In Xi I = kj we have that for each i < d 



Thus, for each i 



I In(xj) — ln(xj)| > I lne| — ki. 
,d, 

^ m 

AxA < — — ; — > c,-|Ax,-|. 



j=d+l 



L&tAmax = snpj^^a+i,...,m}^\^^j\} and 5(e) = sup,g|i^ ,i}(|lne| - h) \ We 
know Amax 7^ because if it were equal to zero, then the right hand side of 
equation (fTTI) would be zero, which it can not be as it is strictly larger than the left 
hand side. Combining the above shows that for each i = 1, . . . , d, 

\Axi\ < 6{e)CA 

max ■ 

Now we consider the vector Ax = x(T^) — x{t^) E S. Normalizing the vector Ax 
by dividing each entry by A^ax then produces a vector vie) = -r-^Ax with the 
following properties: 

1. v{e) E S. 

2. For:i= l,...,d, \vi{e)\ < 5{e)C. 

3. There is at least one entry in t>(e) with norm 1 (the one for which the maxi- 
mum in the definition of Amax was achieved), and none have a higher norm. 
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4. 1 < \v{e)\ < m. 

Property 4 follows from property 3. e > was arbitrary, so we may consider 
a sequence {e„} such that e„ > e„+i and e„ 0. For each e„ we may redo 
the work above. This leads to a sequence of vectors {f (e„)} and a sequence of 
numbers {5(e„)} such that 5(e„) and for each n all four properties above 
hold. Because each vector from the sequence {f (e„)} is contained in the compact 
space {v : \ < \v\ < m} fl S, there is a convergent subsequence (e„^)} and a 
vector zq such that f (e„^) zq & {v : 1 <\v\ < m} n S* C S", as A; ^ oo. Note 
that zq can not be the zero vector because |2;o| > 1. However, 5{tn^) 0, and so 
the first d components of zq are equal to zero. Hence, supp(^o) C supp(?/). □ 

Theorem 3.3. Let {S, C, 71} be a weakly reversible, deficiency zero, mass action 
chemical reaction system with dynamics given by equation ([T]). Suppose that for 
each c G M>q and each semi-locking set W, the set [(c + 5*) fl M>q] fl Lw is either 
empty or discrete. Then the unique positive equilibrium of each stoichiometric 
compatibility class guaranteed to exist by the Deficiency Zero Theorem is globally 
asymptotically stable relative to its compatibility class. 

Proof. We suppose, in order to find a contradiction, that there is a positive equi- 
librium, X, that is not globally asymptotically stable relative to its compatibility 
class. By Theorems 12. 5[ 12.71 and |2.9[ there is a semi-locking set W, a ^ G M^g, 
and a vector y such that y E [{x + S) H M>q] fl L^/ and y G By Propo- 

sition |X2l there exists a nonzero zq E S such that supp(2;o) C supp(y). Because 
y E X + S and zq G 5, for any > we have y + rfZQ E x + S. Further, because 
supp(zo) C supp(y), if 7] is small enough we have that y + r]ZQ G M>o fl Lw- But 
this is valid for all r] small enough, and so [(s + S*) fl M>q] fl Lw is not discrete. 
This is a contradiction and so the result is shown. □ 

Corollary 3.4. Suppose that for a weakly reversible deficiency zero chemical re- 
action system with mass action kinetics, each semi-locking set is a locking set. 
Suppose further that within each stoichiometric compatibility class, the set of 
equilibria on the boundary is discrete. Then the unique positive equilibrium of 
each stoichiometric compatibility class guaranteed to exist by the Deficiency Zero 
Theorem is globally asymptotically stable relative to its compatibility class. 

Proof. Because each semi-locking set is a locking set, the set of boundary equilib- 
ria for a given compatibility class is precisely given by Uvk[('^ + 5') fl ]R>o] ^ ^w-, 
where the union is over the set of semi-locking sets. Therefore, each [(c + 5*) fl 
M>q] n Lw is discrete and invoking Theorem 13 . 3 1 completes the proof. □ 
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Corollary 3.5. Suppose that a weakly reversible deficiency zero chemical reaction 
system with mass action kinetics has only one linkage class. Suppose further 
that within each stoichiometric compatibility class, the set of equilibria on the 
boundary is discrete. Then the unique positive equilibrium of each stoichiometric 
compatibility class guaranteed to exist by the Deficiency Zero Theorem is globally 
asymptotically stable relative to its compatibility class. 

Proof. For single linkage class systems that are weakly reversible, each semi- 
locking set is a locking set. Using Corollary 13 .41 completes the proof. □ 

3.1 Connection with extreme points 

We connect our results to a condition on the extreme points of the positive stoi- 
chiometric compatibility classes. 

Proposition 3.6. For y G R>o, let W = {Xj : ?/j = 0} = supp{y)'^ . Then the 
following are equivalent: 

( i) y is an extreme point of (y + 5*) fl M>q. 

(ii) [{y + S)^Wl,]i^Lw = {y}. 

Proof. ((0) =^ dn])) Suppose that (0) is true, but that du]) is not. Then, because 
\{y + S") M>o] n is not discrete, there exists a G S* Ly^ such that for 
sufficiently small t,y±tv e [{y + 5) H K>o] fl ^w- Noting that y = {l/2){y + 
ef ) + (1/2) (?/ — ev) then shows y is not an extreme point of {y + S)r\ M>q, which 
is a contradiction. 

(® =^ ©) Suppose that dn]) is true, but that (Q) is not. Because y is not an 
extreme point of {y + S) D M>0' there exists nonzero vectors vi ^ y, V2 ^ y in 
{y + S)f] M^o and < A < Tsuch that 

y = \vi + il- \)V2. (12) 

Because vi, V2 E (y + S) f]]R>Q, there exists u,w E S such that vi = y + u 
and f 2 = y + w, and Ui, > if Xj G W. However, because A, 1 — A > 0, 
we see by equation (fT2l) that Ui,Wi = for all Xj G W. Therefore, vi,V2 G 
[{y + S)f] W^o] n Lw, contradicting ©. □ 

Theorem [33] can now be reformulated in the following way. 
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Theorem 3.7. Let {S, C, 71} be a weakly reversible, deficiency zero, mass action 
chemical reaction system with dynamics given by equation ([T]). For any boundary 
pointy G [(c+S*) n]R>o] of a positive stoichiometric compatibility class, letWy = 
{Xi : Hi = 0} = suppiy)'-" . Finally, suppose that Wy is a semi-locking set only if 
y is an extreme point. Then the unique positive equilibrium of each stoichiometric 
compatibility class guaranteed to exist by the Deficiency Zero Theorem is globally 
asymptotically stable relative to its compatibility class. 

Proof. Suppose 14^ is a semi-locking set. Let y E be such that there exists a 
c G M™o with y ec + S. If no such y and c exist, then [(c + S)f] M^q] f]Lw = ^ 
for all c G M™q. If such a y and c do exist, then W = Wy and, by assump- 
tion, y is an extreme point. Thus, by Proposition 13. 6[ [{y + S) n^>o] fl = 
[{c + S)(~] M>o] n is discrete. Invoking Theorem 13.31 completes the proof. 

□ 

4 Examples 

We begin with an example found in [6] for a receptor-ligand model. See (SI for 
full details. 

Example 4.1. Consider the following system, which we assume has mass action 
kinetics 

2A + C ^ A + D 

n n ■ (13) 

B + C ^ E 

For this example there are four complexes, one linkage class, and the dimension 
of the stoichiometric compatibility class is easily verified to be three. Therefore, 
the system has a deficiency of zero and our results apply. The minimal semi- 
locking sets (that is, those that must be contained in all others) are given by Wi = 
{A, B, E}, W2 = {A, C, E}, and W3 = {C, D, E}. Therefore, showing that the 

Lw, is discrete for any c G M>o would also show that the 
sum over all semi-locking sets is discrete. For this example, it is easily verified 
that 



S = Span < 
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One method to show that for a given i the set [(c + S") fl M>o] ^ is at most 
discrete is to demonstrate that there are no non-zero vectors contained in S with 
support given by Wf. This method bypasses the need to check whether the in- 
tersection [(c -|- 5*) n M>q] n is non-empty. It is easily verified that there are 
no non-zero vectors contained in S with support given by W^' or In [0 it is 
shown that for each c G M>0' (c + S) fl M>o ^^^^ intersect one of Lyy-^ or Lw2- 
Therefore, there are always equilibria on the boundary, however, by our results or 
those found in [|6J, they will not affect the global stability of the interior equilibria. 

Let t/3 = {x G I supp(x) G W^]. It is easy to show that U^{^S = 
span{[2, —1, 0, 0, 0]^}. Thus, the method used in the previous paragraph does not 
work. Therefore, for our results to apply, we need to verify that [(c + S') fl M>q] fl 
Lvi/3 = 0, for any c G M5,o- Because is characterized by having the last three 
entries equal to zero, in order to prove that [(c + 5) fl M>q] fl = for any 
it is sufficient to show that the space spanned by the last three entries 



c G 



of the vectors in (|T4l) does not contain a vector with strictly negative components. 
We have 
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(15) 



which does not include a strictly negative vector. Thus, [(c + 5) fl ]R>q] fl = 
for any c G M>o- Combining all of the above with Theorem 13.31 shows that for 
any choice of rate constants and initial condition, the system (fT3] ) has a globally 
asymptotically stable equilibrium value. 

Example 4.2. Consider the system 

2A ^ A + B ^ B + C. (16) 

There are three complexes, one linkage class, and the dimension of the stoichio- 
metric compatibility class is two. Therefore, the system (fT6l) has a deficiency of 
zero. The minimal semi-locking sets are Wi = {A, B} and W2 = {A,C}. The 
stoichiometric subspace is of dimension two and the quantity A + B + C is con- 
served. Thus, each stoichiometric compatibility class is a plane that intersects 
each of Lvki = {"^ ^ "^i = "^^2 = 0, 7^ 0} and = {v ■ Vi = = 0,V2 0} 
in precisely one point. See Figure |4T1 Therefore, by Theorem l3.3[ for any choice 
of rate constants and initial condition, the system (fT6l) has a globally asymptoti- 
cally stable equilibrium value. We note that it is easily verified that the eigenvalues 
of the linearized problem around the equilibria associated with the semi-locking 
set Wi are all zero, and so the results of flU do not apply here. 
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Figure 4.1: The stoichiometric compatibility class for the system (fT6l) is a plane 
defined by the conservation law A + B + C = M for some M > 0. This plane 
intersects the sets L^-^ and Li^^ (bolded axes) each in precisely one point. 

Example 4.3. Consider the system 

2A ^ A + B 2B ^ A + C. (17) 

There are four complexes, two linkage classes, and the dimension of the stoichio- 
metric compatibility class is two. Therefore, the system (fTTI) has a deficiency 
of zero. The only minimal semi-locking set is W = {A, B}, and this is also a 
locking set. The stoichiometric subspace is of dimension two and the quantity 
A + B + C is conserved. Thus, each stoichiometric compatibility class is a plane 
that intersects Lw = {v : vi = V2 = 0,V3 ^ 0} in precisely one point. There- 
fore by Theorem l3.3l or Corollary 13. 4[ for any choice of rate constants and initial 
condition, the system (fTTI) has a globally asymptotically stable equilibrium value. 
It is easily verified that the boundary equilibria are not hyperbolic with respect to 
their compatibility classes, and so the results of [6} do not apply in this case. 

5 Non-mass action kinetics 

In [fT4l . Sontag extended the Deficiency Zero Theorem to systems with non-mass 
action kinetics. He considered weakly reversible, deficiency zero system whose 
kinetic functions are given by 

where the functions 6'j : M — > [0, oo) satisfy the following: 
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1 . Each 6i is locally Lipschitz. 



2. 9i{0) = 0. 

3. I \n{ei{y)\dy < oo. 

4. The restriction of 9i to M>o is strictly increasing and onto M>o. 

To prove the local asymptotic stability of the unique equilibrium within each sto- 
ichiometric compatibility class the following Lyapunov function was used 



where Pi{s) = In^^j(s) and x is the unique equilibrium within the positive stoi- 
chiometric compatibility class. Note that 9{x) = \x\ gives mass action kinetics, 
in which case the Lyapunov function given in equation (fT9l ) is the same as that in 
equation The only dynamical property of the Deficiency Zero Theorem used 
in this paper was that W{x) — oo as Xi 0. We note that for the Lyapunov 
function (fT9l) . we still have that property because 



and pi{xi) = \n9i{xi) — oo as Xj — by the properties of 9i{-) given above. 
Therefore, our results in this paper, and in particular Theorem l3.3l Corollary 13. 4[ 
Corollary 13. 5[ and Theorem 13 .7 [ are valid in the setting (fTSi) . 
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